//report fit of model with "other" moments, like LF frictions and lovenheim/hilger moments
import delimited "$dir/Model/simulated_data/simulated_data_base.csv", clear

//renaming
ren v1 theta
ren v2 e
ren v3 x
ren v4 a
ren v5 t
ren v6 p
ren v7 I
ren v8 H 
ren v9 tau
ren v10 s
ren v11 s_val
ren v12 prob_diff_I
ren v13 prob_diff_H
ren v14 l 
ren v15 m
ren v16 w
ren v17 tax
ren v18 c
ren v19 util
ren v20 switch
ren v21 gov_exp
ren v22 sim
ren v23 i_v
ren v24 debt
ren v25 payment

su debt if t == 25
su debt if t == 25 & debt>0
tab e if t == 25

gen debt_present = (debt>0)

su debt_present if t == 25
su debt_present if t == 25 & i_v<4
su debt_present if t == 25 & i_v>=5 & i_v<10
su debt_present if t == 25 & i_v>=10 & i_v<15
su debt_present if t == 25 & i_v>=15


tab debt_present if t == 30 
tab debt_present if t == 32 & e == 3
tab debt_present if t == 33 & e == 3
tab debt_present if t == 34 & e == 3
tab debt_present if t == 35
tab debt_present if t == 40
tab debt_present if t == 45
tab debt_present if t == 50




//report fit of model with "other" moments, like LF frictions and lovenheim/hilger moments
import delimited "$dir/Model/simulated_data/simulated_data_base_nolim.csv", clear

//renaming
ren v1 theta
ren v2 e
ren v3 x
ren v4 a
ren v5 t
ren v6 p
ren v7 I
ren v8 H 
ren v9 tau
ren v10 s
ren v11 s_val
ren v12 prob_diff_I
ren v13 prob_diff_H
ren v14 l 
ren v15 m
ren v16 w
ren v17 tax
ren v18 c
ren v19 util
ren v20 switch
ren v21 gov_exp
ren v22 sim
ren v23 i_v
ren v24 debt
ren v25 payment

su debt if t == 25
su debt if t == 25 & debt>0
tab e if t == 25
